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ABSTRACT 

The search for the onset of star formation in pre-stellar cores has focussed on the 
identification of an infall signature in the molecular line profiles of tracer species. 
The classic infall signature is a double peaked line profile with an asymmetry in the 
strength of the peaks such that the blue peak is stronger. L1689B is a pre-stellar core 
and infall candidate but new JCMT HCO + line profile data, presented here, confirms 
that both blue and red asymmetric line profiles are present in this source. Moreover, a 
dividing line can be drawn between the locations where each type of profile is found. It 
is argued that it is unlikely that the line profiles can be interpreted with simple models 
of infall or outflow and that rotation of the inner regions is the most likely explanation. 
A rotational model is developed in detail with a new 3D molecular line transport code 
and it is found that the best type of model is one in which the rotational velocity 
profile is in between solid body and Keplerian. It is firstly shown that red and blue 
asymmetric line profiles can be generated with a rotation model entirely in the absence 
of any infall motion. The model is then quantitively compared with the JCMT data 
and an iteration over a range of parameters is performed to minmize the difference 
between the data and model. The results indicate that rotation can dominate the line 
profile shape even before the onset of infall. 

Key words: radiative transfer - submillimetre - stars: formation - ISM: clouds - ISM: 
molecules - ISM: individual: LI 689B 



1 INTRODUCTION 

The physics of the collapse process that leads to the for- 
mation of a star is still not understood. Although many 
models of star formation have been proposed, observational 
difficulties have limited the degree to which they can be 
tested. Molecular line profiles, emitted in the mm and sub- 
mm regimes, potentially offer the best opportunity to ex- 
tract dynamical information about the collapse process that 
leads to the formation of stars. For some species the abun- 
dances are great enough that self- absorption of the line takes 
place and if the gas is undergoing infall then the line can as- 
sume a ch aracteristic d ouble-peaked profile with a stronger 
blue wing. lEvansI 1^)99) has reviewed how such line profiles 
are generated and this is discussed further below. Much ob- 
servational effort has been concentrated on detecting and 
characterising blue asymmetric line profiles in the hope of 
using them to extrac t dy namical informatio n. For exam- 
and lZhou et alJ ll_993j ) mod elled the 
dl977t) 'inside- 



pie, IChoi et al.1 <|1995 

pre-stellar core B335 as undergoing a Shu 
out' collapse. However, the observational data in these cases 
consisted of one or a few molecular lines at a single po- 
sition in the cloud. Though a consistent fit to the data 



could be made using a Shu collapse model there was no 
guarantee that the fit was unique and lWilner et all (l2000t) 
showed that such a model was not able to explain their 
CS (5 — 4) data of B335. While the search for infall in pre- 
stellar cores has become focussed on identifying blue asym- 
metric line profiles, it is becoming clear that the interpreta- 
tion of the molecular line profiles from collapse candidates is 
much more complex than previously thought with the shape 
of the line profile being highly sens itive to the chemical and 
physical properties of the gas (e.g . iRawlines fc Yate'sll200ll : 
IWard-Thompson fc Bucklevll200l|i . 

In this paper, an additional possible cause of asymmet- 
ric line profiles is considered: rotation, which is suggested 
by the observed spatial separation of red and blue asym- 
metric profiles in JCMT line profiles obtained from the pre- 
stellar core L1689B, reported below. L1689B is modelled 
as a system in which the central regions are dynamically 
active and undergoing rotation while the outer regions are 
quiescent. A 3D molecular line transport code is employed 
to model the core and to directly compare the results with 
recently obtained JCMT data. In Section 2 observations of 
L1689B are reviewed and the JCMT data that is to be mod- 
elled is presented. The model and code are described in Sec- 
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tion 3. The results are presented and discussed in Section 4 
an d compared wi t h oth er recent studies, particularly that 
of iBelloche et al.l (|2003) who have recently demonstrated 
the presence of rotation in the more evolved Class source, 
IRAM 04191+1522. It is concluded in Section 5 that the 
asymmetric line profiles can be very well modelled if L1689B 
is undergoing rotation and that therefore rotation dominates 
any infall motions in this source. 



2 THE PRE-STELLAR CORE L1689B 

L1689B is a pre -stellar core located in Ophiuchus. 
IShirlev et all feOOOj) present SCUBA maps of L1689B that 
confirm that there is no pro tostellar object in the centre of 
the core. lEvans et al] i200 1]) use a self-consistent dust emis- 
sion model to constrain the properties of L1689B and con- 
clude that the dust temperature decrea ses from the edge to 
the centre, a conclusion also rea c hed b v lAndre et al] dl996f) 
and IJessop fc Ward-Thompson! i200l]) . The density struc- 
ture of the core is mildly flattened in the centre and de- 



lEvans et al] J200ll 


using a Bonner-Ebert sphere density dis- 


tribution (see e.g. 


McLauehlin & Pudritdll996l lAlves et al 


2001 


: Evans et a 


J 120011: IWhitworth & Ward-ThomDSon 


2001 


; Wuchtcrl & Tscharnuter 200l| for recent work involv- 



ing Bonncr-Eb ert spheres ). 

iRedman et'a l. (2002) carried out C 17 J = 2 -> 1 
observations towards the pre-stellar core L1689B. By ex- 
amining the relative strengths of the hyperfine components 
of this line the optical depth was calculated. This allowed 
accurate CO column densities to be determined. The hy- 
drogen column densities that these measurements imply are 
substantially smaller than those calculat ed fr om SCUBA 
dust e mission data of lShirlev et al] fcOOOl) and lEvans et alJ 
<200ll) . Furthermore, the C 17 J = 2 — » 1 column den- 
sities are approximately constant across L1689B whereas 
the SCUBA column densities are peaked towards the cen- 
tre. The most likely explanation is that CO is depleted 
from the central regions o f L168 9B, as also suggested by 
IJessop fc Ward- Thompson! i200l]) . Evidence of CO deple- 
tion has also been found in several othe r prestellar cores 
JCaselli et alj[l 999: Bacmann et al.l l2002t |Jorgensen^^l" 
l2002l:|Tafalla et al.l2002l) . It was estimated bv lRedman et al 



( 2002) that within about 5000 AU of the centre of L1689B, 
over 90% of the CO has frozen onto grains. This level of 
depletion can only be achieved after a duration of time that 
is at least comparable to the free-fall timescale. While the 
CO emitt ing gas appears to be quie scent, line profiles ob- 
tained by iGregersen fc Evansl ]200Cf) show the presence of 
both blue and red asymmetric line profiles in this source. 
New JCMT observat ions, described below, confi rm the ob- 
servational results of IGregersen fc Evansl (12000) and reveal 
an underlying order to the locations of the blue and red 
asymmetric profiles. 



3 JCMT OBSERVATIONS OF L1689B 

The bulk of the observations were carried out using the 
heterodyne receiver RxA3 at the James Clark Maxwell 
Telescope (JCMT), Mauna Kea, Hawaii on the nights of 



2001 Feburary 24-26, 2001 August 19-21, 2003 May 3-5. 
Additional data were collected on various dates in 2001. 
HCO+ J = 3 -» 2 (267.5665 GHZ) line profiles were ob- 
tained for seventeen positions across the centre of L1689B. 
Some of the data points were collected as part of a larger 
survey of infall candidate objects, the results of which will 
be reported in Redman et al. (2004) . The data were reduced 
in the standard manner using the SPECX software package 
and were converted to the T m b scale using a main beam 
efficiency of 0.7. 

The seventeen positions ar e at different offsets and a t 
shorter spacings than those of IGregersen fc Evansl {2000). 
The line profiles are displayed in FigQ Again, the presence 
of both blue and red asymmetric line profiles is apparent. 
The two datasets cover different, though overlapping, parts 
of the cloud (note that, as is often the case for prestellar 
cores, the position of the centre of the cloud is somewhat un- 
certain - see, e.g., the SCUBA dust emission map of Shirley 
et al 2000). While in our map the majority of the profiles are 
red asymmetric line profil es, the opposite is the case for the 
IGregersen fc Evans! JgOOO jl dataset. Where the pointing po- 
sitions between the IGregersen fc Evansl <l2000l) dataset and 
the present study overlap (in the north-east [top right] cor- 
ner) it was verified that the datasets match each other . 

Inspection of the data of lGregersen fc Evansl i2000f) and 
Fig reveals that the blue and red asymmetric line profiles 
are found on either side of an axis that runs from SW to NE. 
This axis is displaced to the NW of the zero offset position 
of our dataset but would run through the centre of the com- 
bined datasets. The change in asymmetry is in br oad accord 
with t he lower spatial resolution CS 2 — ► 1 data of lLee et al] 
( 1999) (obtained with FCRAO) in which the sense of asym- 
metry of the line profiles switches from blue in the north of 
L1689B to red in the south. This indicates that there is dy- 
namical activity in the HCO + emitting gas. Since the outer 
layers of the cloud, as traced by C 17 Q, are quiescent (the 
hyperfine structure is well resolved - IRedman et al] 2002) 
the dynamically active gas must lie within the CO depleted 
region. 



4 MODELLING 

A blue asymmetric line profile in an infalling core will be 
formed if (i) the infall velocity increases less slowly than 
1/r and (ii) the excitation t emperature increases towards 
the centre (e.g. lEvanslll999 | ). The dust rad iative transfer 
modelling carried out bv lEvans et al] (1200 it) demonstrates 
that the second of these conditions may not be met in an 
object like L1689B since the dust temperature (and, if they 
are well coupled, the gas temperature) decreases towards 
the centre of the cloud. Furthermore, of course, the presence 
of red asymmetric line profiles as well as blue asymmetric 
profiles toward the centre of the cloud means that simple 
infall cannot explain the observational data. 

Since star formation is usually accompanied by out- 
flows, one possible explanation is that an outflow is the cause 
of the changes in the sense of the asymmetry. However, this 
can be discounted here since the CO would be kinemati- 
cally disturbed and the sharply resolved hyperfine structure 
would be blended away. Furthermore, in a starless core with 
an interior depleted in heavy molecules, such as L1689B, an 
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Figure 1. HCO + J = 3 — > 2 line profiles from positions across the centre of L1689B. The coordinates of the pointing centre are RA 16 
31 47.0, DEC -24 31 45.0 (B1950). The axis dividing blue and red asymmetric profiles is displaced to the NW of the zero offset position 
of our dataset 



outflow is not yet expected to have been initiated. Instead, 
we model the system as consisting of a quiescent outer halo 
where the CO is emitted, with a rotating inner HCO + emit- 
ting region. Figure is a sketch of the model. 



4.1 General description of the code 

In dark molecular clouds, the density is usually too low for 
local thermodynamic equilibrium (LTE) to apply, the opac- 
ity is usually too high for an optically thin approximation 
and the systemic velocities are too low for the large veloc- 
ity gradient (LVG) or Sobolev approximations to be valid. 
We therefore generate model spectra and intensity maps us- 
in g a non-LTE num erical radiative tr ansfer code (described 
in iKeto et alj|2003l and also used by iRawlines et~alll200-tfl 
The code is fully 3 dimensional and em ploys the accelerated 
lambd a iteration (ALI) algorithm of iRvbicki fc Hummer! 
Jl99lT) to solve the molecular line transport problem within 
an object of arbitrary (three-dimensional) geometry. This 
repre sents a cons iderable improvement over the early mod- 
els of lKetol dl99CD . 

The local line profile is specified by systemic line-of- 
sight motions together with thermal and turbulent broad- 
ening. In addition to essential molecular data (such as colli- 
sional excitation rates) the physical input to the model con- 




Figure 2. Sketch of the model of L1689B. Th e radius of the 
rotating centre is 3000 AU. Re dman et alj |2002|> estimated that 
most of the CO is frozen onto grains within the central 5000 AU 
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sists of the three-dimensional velocity, density and temper- 
ature structures. The radiative transfer equations are then 
numerically integrated over a grid of points representing sky 
positions. The main source of uncertainty in the radiative 
transfer calculations are the collisional rate coefficients, but 
this should only be manifest in the uncertainties in the line 
strengths (which may be inaccurate by as much as 25%); 
the line profile shapes and and relative strengths should be 
less affected. 

For our purposes we define a spherical cloud within a 
regularly spaced cartesian grid of 30 x 30 x 30 cells. Within 
each cell, the temperature, H2 density, molecular species 
abundance, turbulent velocity width, and bulk velocity are 
specified. Externally, the ambient radiation field is taken to 
be the cosmic microwave background. At a given viewing 
angle to the grid, the line profile is calculated at specified 
offsets from the centre by integrating the emission along 
the lines of sight. It is then possible to regrid the data in 
the map plane and convolve with the (typically Gaussian) 
beam pattern of the telescope, which for the JCMT is 15". 
The spectra are computed with a frequency resolution of 
0.01 km s _1 . 



4.2 Comparison with observational data 

Once a plausible general model has been set up, observa- 
tional data can be read directly into the code and multiple 
models with varying parameters can be calculated and the 
difference between the model and the data minimized. The 
dimension of the model cube is selected so that the data 
points are positioned at integer values of the cube spacing. 
The data is regridded so that the data channels match those 
of the model output. After generating a model cube, the 
code compares the model with the data and calculates \ 2 
for that fit. The numerical techniques of simplex and simu- 
lated annealing are used to drive the refining process where 
the code adopts new pa rameters and calc ulates a new cube; 
full details are given in lKeto et al.l i2003fl . A typical fitting 
run may require the calculation of several hundred model 
cubes. 



4.3 L1689B model 

The model parameters are listed in Table 14*31 along with ei- 
ther the fixed adopted values used in the code or the range 
over which a parameter was allowed to vary. The number 
of model runs required to explore the parameter space in- 
creases rapidly with the number of free parameters so those 
parameters which have the best empirical contraints were 
given fixed values. 

The density profile used is that of a Plummer sphere 
which is defined as 



Parameter 



Adopted value or range (best fit value) 



p(r) 



■ 2 + rl 



(1) 



where po is the central density and ro is the Plummer ra- 
dius. This simple density law has the property that as r — > 0, 
p(r) — > po while for large r, p(r) — ■> po/r 2 . The outer radius 
of the cloud is defined as 6000 AU. The Plummer profile is 
also a reasonable approximation to (and easier to deal with 
than) the Bonner-Ebert sphere density distribution. This 
latter distribution appears to describe the density profile of 



Outer radius 
Plummer radius ro 
Plummer density po 
Temperature 
HCO+ abundance 
Projected velocity 
Rotation radius 
Turbulent velocity 



6000 AU 
1000 AU 
10 6 cm" 3 

10 K (edge) - 7 K (centre) 

5.0 x lO" 10 - 1.0 x 10~ 8 (7.84 x 10~ 9 ) 

0.1 - 0.4 km s- 1 (0.121) 

3000 AU 

0.12 - 0.18 km s" 1 (0.157) 



Table 1. Model parameters used. For the parameters that were 
allowed to vary, the best fit value is given in brackets 



pre-stellar cores and is begining to be used routinely in mod- 
elling work. The adopted values of ro and po were 1000 AU 
and 10 6 cm -3 respectively. These values were chosen to ap- 
proximate the Bonner-Ebert sph ere fit to the density distri- 
bution of L1689B carried out bv lEvans et alJ (|200 if) . Simi- 
larly, an em pirical approximati on to the dust temperature 
profile fit of lEvans et al.l i200ll) was employed to trace the 
gas temperature (the gas density is sufficiently high that the 
two temperatures can be considered well coupled). The tem- 
perature is 10 K in most of the cloud and then drops, over a 
short distance, to 7 K in the centre. In fact, trial runs with 
a constant temperature of 10 K throughout the cloud yield 
very similar results to this temperature profile. 

The fractional abundance of HCO + was taken to have 
a constant value from the range 5.0 x 10 -10 — 1.0 x 10 -8 . 
The assumption of a constant abundance is of course a crude 
approximation and a self-consistent chemistry should be in- 
cluded in the cl oud. However, the ch emistry of HCO + is 
complicated and lRawlings et alJ dl992fl show that its abun- 
dance intially rises as depletion ensues. The sticking coef- 
ficient of HCO + on to charged and neutral dust grains is 
also not well known. Therefore an attempt to include the 
chemistry does not yet seem worthwhile. 

A simple solid body rotation law was adopted for the in- 
ner regions of the model cloud. A rotation axis that matches 
the dividing line between the red and blue asymmetric pro- 
files was imposed. The range of projected velocities inves- 
tigated is 0.1 — 0.4 km s - . The rotation radius is an es- 
timate of where the inner regions stop rotating and is es- 
timated as being at 3000 AU (c f. 5000 AU for the CO de- 
pleted 'hole' lRedman et al .120021) . Finally, the uncertain tur- 
bulent velocity width was allowed to vary over the range 
0.12 - 0.18 km s _1 . 



5 RESULTS AND DISCUSSION 

The full set of synthesised HCO + J = 3 — » 2 line profiles 
for the best fit model is displayed in Figure [3] For clarity, 
the rotation axis points vertically. The line profiles show 
clear asymmetries with the sense of the asymmetry switch- 
ing between red and blue on either side of the rotation axis. 
The self-absorption is strongest towards the centre of the 
model cloud and the asymmetry is strongest there. Away 
from the centre of the cloud, the drop in column density 
and self-absorption results in a reduced line strength and 
less pronounced asymmetry (the difference in strength be- 
tween red and blue peak; the relative magnitude of the ab- 
sorption trough). These results confirm that it is possible to 



© 2004 RAS, MNRAS 000,[T]|7| 



Rotation of L1689B 5 



Red asymmetric 



Blue asymmetric 



A 



JX. 



JX. 



A 







Jl 



JA 



JA 



JA 



XX 



Jl 



A 



Jl 



ji 



A 



H 



A 



H 



JL 



Jl 



XX. 



Jl 



al 



.XL 



Jl 



JA 



-XL 



Jl 



JA 



_Q_ 



Jl 



A 



al 



ji 



A 



-1.S0.9D.DQ.5U 



Jl. 



A 







Jl 



Jl 



XX 



A 



ji 



A 



al 



A 



ji 



A 



ji 



A 



JL 



-0_ 



A 



Jl 



xx 



Jl 



Figure 3. HCO + J = 3-^2 line profiles from the rotation model. For clarity the array of profiles is displayed with the rotation axis as 
shown. The line profiles switch between red asymmetric and blue asymmetric on either side of the rotation axis. 



generate both red and blue asymmetric line profiles entirely 
in the absence of any infall or outflow motions. 

Figure 0] is an overlay of the best fit model described 
above with the JCMT data of Redman et al (2004). The 
rotation axis is marked on the figure. Given the simplicity 
of the model, the fit is excellent with the sense of the asym- 
metry in the line profiles, their shape and the fall-off in line 
strength away from the centre being reproduced well. The 
most obvious of the discrepancies is that the line strengths 
are overestimated for the brightest lines near the centre of 
the cloud. One explanation for this may be due to the freez- 



ing out of HCO + onto grains. Any freeze-out of the HCO + 
is likely to occur first in the centre of the cloud, where the 
density is highest, and progress outwards. If there is signif- 
icant freeze-out in the very centre of the cloud, this would 
reduce the strength of the line profiles for small impact pa- 
rameters. As noted above, the chemical behaviour of HCO + 
in a core undergoing depletion is complicated. 

iBelloche et alJ i2003f) have carried out a detailed study 
of the Class protostar IRAM 04191+1522 and conclude 
that the inner region of IRAM 04191+1522 is undergoing 
fast differential rotation while the outer regions are only 
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Figure 4. HCO + J = 3 — > 2 line profiles from L1689B. The spacing between the cells is 10". North is up and east is to the right. Overlaid 
are synthesised line profiles generated by the rotation model for the source. 



slowly rotating. Since they also argue that this source has 
only recently become a protostar it is constructive to com- 
pare it with L1689B which may be on the verge of col- 
lapse. Firstly, under the assumption that both objects have 
been adequately modelled, the angular velocity is compara- 
ble: lot = 9.4 km s" 1 at 3000 AU in the best fit model for 
L1689B while for IRAM 04191+1522, ujr = 9±3 km s" 1 at 
2800 AU. IRAM 04191+1522 may be decou pling from the 
outer envelope at a radius of 2000-4000 AU llBelloche et alJ 
2003) which compares with the rotation radius of 3000 AU 
estimated here for L1689B. Solid body rotation is firmly 
ruled out for the inner region of IRAM 04191+1522 with 
differential rotation providing a much better fit to the obser- 
vations (specifically, the angular velocity gradient is found 
to change from approximately lu oc r~ 2,5 to u oc r _1 in 
the inner regions with very slow rotation in the outer re- 
gions). The model results presented here for L1689B also 
supports an u> oc r _1 angular velocity profile. Such a veloc- 
ity gradient indicates that angular momentum is conserved 
and that magnetic brakin g does not dominate the dynam- 
ics; iBelloche et alJ l)2003|) provide an extensive discussion 
of their results in terms of various collapse models. Thus, 
L1689B exhibits strong similarities with IRAM 04191+1522 
in terms of the rotational behaviour. The principal kinematic 
difference is that IRAM 04191+1522 shows clear evidence of 
infall (and outflow) in addition to rotation yet the signature 



of any such motion in L1689B is completely dominated by 
the rotation. It seems reasonable to suggest that L1689B and 
IRAM 04191+1522 are similar objects that are at slightly 
different evolutionary stages, either side of the onset of in- 
fall. 



6 CONCLUSIONS 

A model of the pre-stellar collapse candidate L1689B has 
been presented where the inner regions of the cloud are 
undergoing rotation while the outer part of the cloud is 
static and quiescent. The model calculations demonstrate 
that blue and red asymmetric line profiles arise naturally 
from such a system even in the absence of infall or outflow- 
ing gas motions. The possibility that rotation is present in a 
pre-stellar object must therefore be addressed if asymmetric 
line profiles are to be used to argue that a core is undergoing 
infall. 

A best fit model was generated by varying several of 
the free parameters within reasonable ranges. An excellent 
fit to the data was obtained and we would therefore argue 
that our model is a plausible one for L1689B. 

This study along with other recent work demonstrates 
that great care needs to be taken in interpreting asymmet- 
ric line profiles as evidence for infall. Whi le this means that 
earlier claims for detection of infall (e.g. IZhou et a.lJll99oh 
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ought now be re-evaluated it does not mean the search will 
be fruitless. A combination of chemical models, radiative 
transfer codes and high quality observations should soon 
lead to the first firm identification of infall and the long 
awaited detailed examination of the competing star forma- 
tion models. 
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